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The first metastable triplet state of HeH+ was found to be present in ion beam experiments, with its lifetime 
estimated to be between hundreds of milliseconds and thousand of seconds. In this work, we use ab initio 
methods to evaluate the radiative lifetimes of the six vibrational levels of the a '^1]+ of HeH+. The transition 
a — >■ X ^E"*" is spin-forbidden, but acquires intensity through spin-orbit interaction with the singlet 
and triplet 11 states. Large scale CASSCF/MRCI calculations using an adapted basis set were performed to 
determine the potential energy curves of the relevant states of HeH+ as well as the matrix elements of the 
dipole and spin-orbit operators. The wave functions and energies of the vibrational levels of the a and 
X states are obtained using a B-spline method and compared to previous works. We find that radiative 
^ I lifetime of the vibrational levels increases strongly with v, the lifetime of the v = state being of 150 s. We 

also analyze the contributions from discrete and continuum parts of the spectrum. With such a long lifetime, 
the a state could have astrophysical implications. 



o 

(N 



=3 

1—5 



o 
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I. INTRODUCTION 



Ph. 

a 

2 ' The helium hydride ion HeH+ is one of the most elementary molecular ions and the first to form in the early 
Ci , universe^. Due to its high relative abundance, HeH+ was predicted to be observable in astrophysical objects such as 
c/3 ' planetary^ and gaseous nebulae, or in metal-poor stars^.. However, no infrared emission from HeH+ molecular ion has 

yet been detected from these objects'*'^, although it has been observed in laboratory plasmas for many yearsfi. The 
c/2 formation of HeH+ is mainly due to the radiative association between He and H+ or between Hc+ and Hi. While 
^|]~» it was always supposed that HeH+ formed in its ground X ^T,'^ state, one should also consider the possible role of 
1^ the first metastable triplet state, a ■^1]+. This state can indeed be populated and will not decay by collisions if the 

plasma density is low. As its radiative decay to the ground state is spin-forbidden, it is thus expected to have a very 

long lifetime. 

' In studies on the dissociative recombination of HcH+, it was shown that the a '^E^ state is responsible for a part 
' of the cross section, so that it must be present in the ion beam. The lifetime of this state was postulated by Yousif 
' et to be longer than the one of He(ls2s '^S), which decays by a relativistic magnetic dipole transition and has 
1^ [ a lifetime of approximately 8000 s. In an experiment on the charge-transfer dissociation of HeH+, Strasser et a& 

■ estimated the lifetime of the triplet state to be in the range of a few hundreds of milliseconds, a much lower value 
probably due to coUisional decay. These two estimations provide lower and upper bounds on the lifetime of this state, 
but the difference is so large that it motivates a theoretical investigation of the lifetime. 

We present in this work the results of large scale ab initio calculations performed with the MOLPRO program 
] suited. These calculations take the spin-orbit coupling between the low-lying singlet and triplet E+ and H states into 
' account, allowing the spin-forbidden a "^2+ — > X ^S"*" dipole transition to occur. This mechanism was already used 
', to estimate the lifetime of the first triplet state of NO"^— li^. 

■ The resolution of the vibrational problem in the considered potential energy curves has been done using a B-spline 
\—{ ' basis set method, which allowed us to estimate the contribution of the continuum states to the lifetime. 

a 

_ _ i 

II. THEORY 
A. Lifetime 

The inverse of the lifetime of an initial excited electronic state \i) is given in the electric dipole approximation 

Tr^ = E^^/-^E^'/KW)P (1) 

where the sum extends over all states |/) with an energy Ef < Ei. The Aij are the so-called Einstein coefficients for 
spontaneous emission and /x is the dipole operator. Eif — Ei — Ef is the energy difference between states \i) and |/). 
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In the Born-Oppenheimer approximation, the total wave function for the initial or the final state is expanded as 
a product of electronic, vibrational and rotational wave functions. These three types of motions are represented by 
the quantum numbers m (electronic), v (vibrational) and J (rotational). In Hund's case (a), the wave function for a 
state \mvJ) is 

"^mvJ = CmASs(i?, r)^mvJAsiR)^Mni^^ '/') (2) 

A is a quantum number associated to L^, the projection on the internuclear z axis of the total electronic angular 
momentum L. E is associated to Sz, the projection on the z axis of the total electronic spin S. The total angular 
momentum isJ = N + L + S, where N is the angular momentum for nuclear rotation. 17 = A + E is a quantum 
number associated with J^. The vibrational equation that the functions ipmvjAsiR) niust satisfy is, in atomic units, 

- l^dji + U„iyJAs{R) - EjnvJASj'^PmvJAsiR) = (3) 



2fi 

where UmvjAs{R) is the electronic energy, UiAsiR), corrected by a centrifugal term originating from the rotational 
Hamiltonian, 

UravJAsiR) = UmAs{R) + ( J(J + 1) ~ fl^ + S {S + 1) - E^) (4) 

and ^ is the reduced mass of the system. Note that there should be an additional (L^ + Ly) term due to the diagonal 
Born-Oppenheimer correction, which we have omitted. 

If we neglect the rotational motion, the dipole transition moment between the initial and the final state is then 
given by 

{iv\ti\fv') = {ivAS\{iAST,\(i\fA'S'E')\fv'A'S') (5) 
while the lifetime of the vibrational level v of the excited electronic state i is given by 

= EE^™./-' ' (6) 



where the Einstein coefficients Aiyjyi are 
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Mv.fv' 



J2EijJ{ivAS\{iASJ:\fi\fA'S'i:')\fv'A'S'}\^ (7) 
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To be exact, the sum over the vibrational levels v' of state / in equation ^ should be understood as a sum if v' is 
a discrete (bound) level or as an integral if v' is a continuum (unbound) level. However, as explained below, we will 
use a discretization method to treat the continuum, so that the sum will run over the discrete levels until convergence 
in equation ([6|). As the ground state is the only state below the a "^E^ state, the sum over / will reduce to only one 
term. 

To calculate the lifetime of state it is thus necessary to know: (i) its transition moments with all the states |/) 
lower in energy, and (ii) the vibrational wave functions (and energies) of states \i) and |/). 



B. Spin-orbit coupling 

The transition a '^E+ — X ^E+ is forbidden at all multipole orders due to the spin selection rule AS* = 0. However, 
it can occur through spin-orbit coupling and the matrix element of interest, 

(X iE+|/x|a3E+)so , (8) 

will be non zero. 

We add to the molecular Hamiltonian H the Breit-Pauli H^'~' perturbation term, given hy^ 

^so ^ « ^ ^'.A • - V E ^ + 2s,) , (9) 
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where i and A denote electrons and nuclei, respectively and a is the fine structure constant. The first term in equation 
is the direct spin-orbit interaction, while the second term is the spin-other orbit interaction. As this perturbation 
mixes the spin and orbital angular momenta of the electrons, the correct quantum number is il = A + E. In this 
representation, the states characterized by A and S split into i7 components, whose symmetry can be determined 
by group theory^^ from the direct product of the spatial and spin symmetry species. The a "^E"*" state will split 
into two components corresponding to the 11 and irreducible representations of the '^^cxiv point group. Using 
standard labeling, these components are denoted $7 = 1 and = 0~, respectively. All the ft components are 
doubly degenerate. However, the diagonal elements of can be shown to be proportional to AE so that the three 
components of the ^E+ state are still degenerate^. 
The selection rules for spin-orbit coupling arei^ 

An = 0; AS' = 0,±1; AA = -AE = 0, ±1 ; E+ < — ^ E" (10) 

In accordance with the last rule, there will be no spin-orbit interaction between the a '^E"'" and the "'^S+ states. 
However, if we take into account higher parts of the spectrum of HeH+, and in particular H states, the transition will 
become possible. 

If we consider a ^H and a ■^H state, we can use the rules pU)) to write the X ^E+ and a ■^E+ wave functions in the 
spin-orbit representation in terms of the unperturbed functions as 

\X iE+ )so=ci|A 1e++)+C2|3Ho+) 

|a3E+ )so =C3|a3E+ )+C4|X-) (11) 
\a 3e+)so = C5|a ^E+) + cel^Hi) + cyl^Hi) 

where the coefficients Ci are obtained by diagonalizing the spin-orbit matrix in the basis of the unperturbed functions. 

The relevant matrix element ([S]) in the spin-orbit representation can be evaluated, provided that the mixing coef- 
ficients and the dipole transition functions are known in the unperturbed basis. Due to the splitting of the triplet 
state into two components, the matrix element ^ is split into two parts, according to the value of fl. However, 
the matrix element {X -'^E^^^l/xja ■^E^Jl) vanishes identically due to the fact that the electric dipole operator cannot 
connect 0"*" with 0~ states. The component a '^E^_ will thus decay through another mechanism such as spin-rotation 
or relativistic magnetic dipole perturbations, which are smaller by several orders of magnitude. We will therefore only 
be able to calculate the lifetime of the 17 = 1 component of a "^E"^. However, as the three components are degenerate, 
the total lifetime of this state will be given by 

r-\a 3e+) = I^Y.EL'\{X '^^Ma '^tv')so\' (12) 

v' 

The matrix element ([5]) can be evaluated for the ft — 1 component using (ITT]) as 

{X iE+ |/x|a 3e+)so = ciCe{X '^+\fi\^n) + cacgf H|/x|a + cacyf Hl/x^H) (13) 

All the dipole matrix element occurring in equation (|13p are non zero for the and fiy components of the dipole 
operator, and the extension of equations ([TT|) and (IT^ to the case of more than one singlet or/and triplet H state is 
straightforward. 

III. RESULTS AND DISCUSSION 
A. Dipole transition function 

We will consider here the lower part of the spectrum of HeH+ , which is composed of states that dissociates either 
into H+ -I- He(lsnZ ^'^L) or into H(riZ) -I- He''"(ls). As the various ^E+ and '^E+ states cannot interact through 
spin-orbit perturbation, it is not necessary to include in the spin-orbit calculations other E states than the X ^E+ 
and the a ^E+. On the other hand, all the singlet and triplet H states will contribute to the matrix element 
We have found (see below) that five ^H and three "^H arc sufficient to describe correctly the transition. The states 
included are given in table U together with their dissociation products, and their adiabatic potential energy curves are 
presented in figure [1] 

All the calculations were done using the MOLPRO program-'^^ and an adapted basis set which consists for each 
atom of the aug-cc-pV5Z (or AV5Z) basis sefl-'lS augmented by [3s, 3p, 2d, If] Gaussian-type orbitals optimized to 
reproduce the spectroscopic orbitals of the He and H excited states (seei2. for details). To obtain the potential energy 
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curves for the electronic states, we performed a state-averaged CASSCF— i2i followed by a CI calculation using an 
active space of five a, ten tt and one 5 orbitals. The Breit-Pauli spin-orbit matrix has been calculated on the basis of 
the unperturbed CASSCF eigenfunctions. To include additional correlation effects, we replaced the CASSCF energies 
(the diagonal elements of the spin-orbit matrix) by the calculated CI energies. This matrix is then diagonalized using 
the state interacting mcthod^^ implemented in MOLPRO. The diagonalization of the spin-orbit matrix corresponding 
to all (A, S) states of table H] provides 32 roots corresponding to 4, 4, 9 and 3 states with VI = 0+,0^,l and 2, 
respectively. 



Symmetry U^as{R = 70) 


Dissociative atomic states 


X 


-2.90324307 


H+ + He(ls2 ^5) 




-2.49996040 


H(ls) + He+(ls) 


1 


-2.12491660 


H(2p) + He+(ls) 


2 


-2.12368473 


H+ + He(ls2p ip°) 


3 


-2.05639837 


^H(3p) + ^H(3d) + He+(ls) 


4 


-2.05616425 


H+ + He(ls3d ^D) 


5 


-2.05456333 


^H(3p) - ^H(3d) + He+(ls) 


1 ^'n 


-2.13282525 


H+ + He(ls2p 3p°) 


2 


-2.12491845 


H(2p) + He+(ls) 


3 


-2.05814410 


H+ + He(ls3p ^P°) 



TABLE I. Molecular states included in the calculations, together with their dissociative products and energy at i? = 70 au. 
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FIG. 1. Potential energy curves of the molecular states included in the calculations. 

The evolution of the dipole transition moment © as a function of the internuclear distance is presented in figure 
m It is of order of 10~^ Cg, in agreement with the fact that the spin-orbit interactions are small in light molecules. It 
also decreases rapidly, corresponding to the fact that the value of the spin-orbit matrix elements under consideration 
vanish in the atomic limit. 

On figure [3] are presented the weight of each of the 11 states taken into account into the dipole matrix element 
{X ^E^^^wj/xja •^I]^t;')so. We see that the matrix element is dominated by the contribution from the first ^11 state. 
The inclusion of the second and third ^11 states modify the matrix element by more than 10 %, but the fourth and 
fifth states bring an additional correction of only 2-3% so that it can be supposed that the inclusion of more 
singlet n states will not affect dramatically the lifetime. The contribution of the first triplet 11 state increases with 
the internuclear distance but does not exceed 10%, while the second and third '^11 states contribute to less than 3%. 
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It is therefore reasonable to assume that the value of the lifetime obtained when considering five and three ^11 
states is correct up to a few percents. 
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FIG. 2. Dipole transition function {X ^E++|/x|a '^T,f)so- 




B. Vibrational analysis of the X and a ^E""" states 

The vibrational analysis of the X and a '^E"'" states was done using a B-spline basis set method^^. 

The ab initio calculation of molecular vibrational spectra involving highly excited vibrational states implies the 
variational resolution of the purely vibrational Schrodinger equation by analytical FBS (Finite Basis Set) or numerical 
DVR (Discrete Variable Representation) approaches. In this work, B-spline basis sets are used as alternative FBS 
method. The flexibility of B-splines has been demonstrated in atomic and molecular calculationa^^iS^ by an accurate 
description of both the bound and the continuum states and their efficiency in the resolution of the nuclear motion of 
molecules has been assessed. In opposition with the DVR, which is constrained by a uniform distribution of the grid 
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points, the B-splines allow more flexibility from the definition of different cavities in which the number of grid points 
can be adjusted. 

The definition of a B-spline basis set on a cavity of size L = a^min — a^max starts with the definition of a sequence 
of N real-valued knots {^i}, satisfying ti < i^+i. Using this knot sequence, a set of N polynomials, all with same 
degree, is defined: the B-splines of order K (and degree K — 1). The recursive definition of the B-splines Bi^k{x) is 
the following: 

u I \ J 1 li ti<x < U+i , . 

Bi,i{x) = < . ^ (14) 

I otherwise 

Bi,k{x) = — ^-^Bi^k^i{x) + J''^'' ^ Bi+i^k-i{x) (15) 

It should be mentioned that we are free to choose the order K, the knot sequence and the interval of the basis set. 
The space that is spanned by a B-spline depends on this degree; the higher the degree, the larger the space that is 
spanned. We used B-splines of order 13 on a grid with Xmhi = 0-5 au and x^ax = 100 a-u divided by 800 equidistant 
knot points. 

Extensive studies of the ground state of HcH+ have been done by Kolos and Peek^^ or Bishop and Cheung^^. We 
reproduce the value for the position of the minimum found by these authors, which is 1.463 au. Kolos and Peek 
find a dissociation energy of 16455.64 cm~^ and Bishop and Cheung obtained a value of 16456.15 cm~^ when Born- 
Oppenheimer diagonal corrections are included. Our result of 16465.5 cm~^ is therefore about 9 cm^^ larger. The 
replacement of the AV5Z by the AV6Z basis set leads to an improvement of less than 2 cm~^. 

The vibrational energies of the ground state of HeH'^ have been studied before using diagonal Born-Oppenheimer 
corrections by Bishop and Cheung^l, or with non-adiabatic and relativistic effects by Stanke et a^. As was found in 
these papers, we obtain 12 vibrational levels. Although we do not reproduce the absolute energies of the levels, the 
energy separation between these levels is correct with an error of less than 3 cm^^, as shown in table |lTl The v = 11 
level, not shown in table|lTl has a binding energy of 1.35 cm~^, to be compared with he value of 1.30 cm~^ found by 
Stanke et al. It should be noted that the use of the AV6Z instead of the AV5Z basis set leads to different absolute 
energies for the vibrational levels, but that the spacing remains the same. 

The a '^E+ state has been studied by Michela^S, and a more accurate study was done by Kolo^ using variational 
wave functions in elliptic coordinates. Both authors found the minimum to be located at 4.47 au, but the dissociation 
energy found by Michels is 661.4 cm~^ while the result of Kolos is 849.0 cm^^. We find an equilibrium position of 
4.452 au, and a dissociation energy of 849.79 cm"'^. 

The only calculations on the vibrational levels of the a '^S]+ state were done using the potential energy curve given 
by Michels completed by an analytical expression at large and it was found that this state supports five vibrational 
levels. The binding energies of the levels v = — 3 can be found in this article and in^i where more precise values 
are presented. We reproduce the energy of the levels up to 1 cm~^. However, using our potential energy curve we 
find that this state supports six vibrational levels, with the v — 5 state being bound by less than 2 cm~^. The total 
energies are given in table Hill together with the binding energies. 

The vibrational energies presented in table IIIII correspond to Hund's case (b) for J = so that the centrifugal 
correction in equation (j?]) vanishes. As we have seen, we will only be able to calculate the lifetime of the J7 = 1 
components of the a state. For these, the centrifugal term also vanishes (see equation (U)), so that the vibrational 
energies presented in table |lll] energies are still valid. 



C. Calculation of the lifetime 



As can be seen from equation ([7]), the Einstein coefficients will depend on the dipole matrix element between the 
initial and final electronic state, and on the overlap of the vibrational functions. We have seen that the dipole is very 
small, but in addition the overlap of the vibrational functions is poor since the two states under consideration have 
their minima located at very different geometries, as illustrated on figure 21 

Using the B-splines, it is a simple task to integrate the overlap between the vibrational functions, multiplied by the 
dipole matrix element. We performed the integration on a grid R £ [1, 100]. For the continuum part of the spectrum, 
it is necessary to sum over all vibrational functions of the pseudo continuum spectrum until convergence. 

On figure O are presented the contributions from the discrete and continuum parts of the vibrational spectrum to 
the Einstein coefficients Aif. It is seen that the convergence is reached with 200 continuum functions. To represent 
this contribution, it is necessary to take into account the fact that the vibrational continuum has been discretized by 
multiplying the Aif by the density of states p{E) = 2/{E{vf+i)—E{vf-i)). For the discrete part of the spectrum, there 
is no density of states, but the Aif should still be multiplied by some value to allow comparison since p{E) depends on 
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V 


Basis 1 


Basis 2 


Stanke et afi^ Bishop and Cheun^^ 





2910.40 


2910.57 


2911.02 


2911.29 


1 


2604.12 


2604.13 


2604.21 


2604.32 


2 


2296.13 


2296.02 


2295.64 




3 


1983.20 


1983.02 


1982.13 




4 


1662.19 


1661.90 


1660.45 




5 


1330.39 


1329.98 


1327.91 




6 


987.71 


987.20 


984.50 




7 


643.06 


642.45 


639.35 




8 


330.72 


330.26 


327.49 




9 


117.51 


117.61 


116.22 




10 


24.92 


24.88 


24.44 





TABLE II. Energy difference Ev — E^+i between two successive vibrational levels of the X ^E"*" state of HeH^ in cm ^ and 
comparison with previous works. Basis 1: AV5Z + adapted basis set; Basis 2: AV6Z + adapted basis set. 



V 


-Etot 


Ehind Ev 







-2.502 969 82 


-664.834 


297.67 


1 


-2.501 613 56 


-367.169 


199.92 


2 


-2.500 702 66 


-167.249 


109.99 


3 


-2.500 201 52 


-57.260 


44.61 


4 


-2.499 998 26 


-12.652 


11.53 


5 


-2.499 945 75 


-1.126 




Unbound 


-2.499 940 62 







TABLE III. Total and binding energies of the bound vibrational levels of the a state of HeH"*", as well as the energy 
difference J5„ — E^+i between two successive vibrational levels. Total energies in atomic units, binding energies in cm~^. 



the energy units we choose. Following^A, we use a density for the bound states given by p{E) ~ l/{E{vf+i) — E{vf)). 
On figure [5l we observe the continuity of the resuhs around the dissociation hmit. 

The hfetime, as weh as the relative contribution from bound and continuum states, are presented in table ITVl We 
observe that the lifetime increases with the vibrational number v. We also see that the contribution of the continuum 
to the lifetime is very small (5%) for Vi — 0, but increases with the value of Vi up to 30% for vt = 5. The contribution 
of the higher lying 11 states will probably reduce the lifetime by a few percent, but our calculations provide an upper 
bound of 150 s on the lifetime of the v = level of the a state in the absence of collisional decay. 

Vi Bound Continuum Lifetime (seconds) 






0.949 


0.051 


149 


1 


0.906 


0.094 


211 


2 


0.786 


0.214 


347 


3 


0.709 


0.291 


738 


4 


0.694 


0.306 


2288 


5 


0.692 


0.308 


14267 



TABLE IV. Relative contributions of the discrete and of the continuum part of the vibrational spectrum to the sum of the 
Einstein coefficients Aif and value of the lifetime of the f2 = ±1 component of a ^E"*" for Vi — — 5. 



IV. CONCLUSIONS 

We have determined the radiative lifetime of the a ^S"*" state of HeH+ using ab initio methods. The decay of this 
state onto the ground X ^E"*" state is spin- forbidden but can occur through spin-orbit coupling. We took into account 
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FIG. 4. Adiabatic potential energy curves of the first triplet a ^S"*" state (above) and of the ground X ^E"*" state (below) of 
HeH+ and position of the bound vibrational levels supported by these states. 
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FIG. 5. Contribution from the discrete (circles) and continuum (crosses) parts of the vibrational spectrum to the Einstein 
coefficients for = 1 as a function of the energy of the final vibrational levels. 



the interaction of the a ^S"*" and X ^T,~^ states with the first five and three states of HeH"'" to estimate the 
dipole transition matrix element {X ^E+|^|a ^S+)so- The vibrational energies and wave functions of the a '^E+ and 
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X ^S"*" states were obtained using a B-spline method and were found to agree well with previous calculations. We 
presented theoretical values of the lifetime of the six vibrational levels of the a '^S"'" state. The lifetime is found to be 
of about 150 s for the d = state and increases rapidly with d, as does the contribution of the continuum states to the 
lifetime. Such a long lifetime suggests that HeH"'" could be present in the a ^S"*" state in astrophysical environments. 
We will investigate the radiative association in this state in a separate work. 
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